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The influence of nearest-neighbor diffusion on the decay of a metastable low-coverage phase (mono- 
■ layer adsorption) in a square lattice-gas model of electrochemical metal deposition is investigated 

by kinetic Monte Carlo simulations. The phase-transformation dynamics are compared to the well- 
established Kolmogorov-Johnson-Mehl-Avrami theory [A. N. Kolmogorov, Bull. Acad. Sci. USSR, 
Phys. Ser. 1, 335 (1937); W. A. Johnson and R. F. Mehl, Trans. Am. Inst. Min. Metall. Eng. 135, 
!>• ; 416 (1939); M. Avrami, J. Chem. Phys. 7, 1103 (1939); 8, 212 (1940); 9, 177 (1941)]. The phase 

transformation is accelerated by diffusion, but remains in accord with the theory for continuous 
nucleation up to moderate diffusion rates. At very high diffusion rates the phase-transformation ki- 
netic shows a crossover to instantaneous nucleation. Then, the probability of medium-sized clusters 
O , is reduced in favor of large clusters. Upon reversal of the supersaturation, the adsorbate desorbs, 

but large clusters still tend to grow during the initial stages of desorption. Calculation of the free 
energy of subcritical clusters by enumeration of lattice animals yields a quasi-equilibrium distribu- 
tion which is in reasonable agreement with the simulation results. This is an improvement relative 
to classical droplet theory, which fails to describe the distributions, since the macroscopic surface 
tension is a bad approximation for small clusters. 
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I. INTRODUCTION 
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Kinetic lattice-gas models are a valuable tool for the modeling of nucleation and growth of adsorbate layersii*^ In 
electrochemical systems, due to the presence of solvents and counter ions, common models are often more severe simpli- 
fications compared to surface science in ultra-high vacuum. It is also more difficult to estimate model parameters from 
first principles. 4,5 Though such estimations might seem desirable, they are not mandatory, and simpler ones, for exam- 
ple by fits to experiments, may be used, at least as a start&i In spite of these complications, successful applications of 
the method to electrochemical systems are encouraging (see, e.g., Rcfs. 6 8). In addition, taking a more fundamental 
point of view, such truly microscopic models - and the equivalent two-dimensional Ising models - may be compared 
to prevalent theories of nucleation and growth like the Kolmogorov-Johnson-Mehl-Avrami (KJMA) theory of phase 
transformationSf2ii2iiiii2ii£ thus enabling one to investigate the latter's validity and underlying assumptions. 14 ' 15,16 
The validity of the KJMA theory has been established not too far below the critical temperature over a wide range of 
magnetic fields (corresponding to supersaturation in the lattice-gas language)i 14 i 15 i 16 However, only spin-flip (adsorp- 
tion/desorption) moves have been included in the dynamics. In the present work, we extend these numerical studies 
to the effects of lateral adsorbate diffusion to nearest-neighbor positions (corresponding to spin-exchange moves in 
magnetic language). 

Our work is also inspired by recent electrochemical pulse experiments on the dynamics of the processes that occur 
after lifting the (22 x \/3) reconstruction of the Au(lll) surfaced In that study, during a short positive pulse of 
the electrode potential, the reconstruction is lifted immediately, and excess Au atoms are released onto the then 
(1 x 1) substrate, where they quickly nucleate and form islands. After the pulse, the surface slowly reconstructs again, 
consuming the excess ions. During this process, however, the adsorbate structure still undergoes post-deposition 
dynamics, where small islands tend to decay faster, and large islands gain in size before they eventually decay. 
However, we emphasize that the connection of our model to these experiments is only qualitative. As the source for 
the adsorbate, we use ions from the solution, whose supersaturation can be pulsed. Moreover, we neglect the geometry 
of the surface and the dynamics of the lifting and reconstruction process. Also, the peak coverages in our simulations 
are much higher than in the experiment. We rather focus on the general characteristics of the growth and subsequent 
decay of submonolayer adsorbates with attractive interactions during and after a pulse of their supersaturation, and 
how these characteristics compare to the KJMA theory. Though we keep the model as simple as possible, nearest- 
neighbor diffusion is an essential characteristic of the experimental system and must be included in the dynamics of 
the model. 

We shall not repeat here the discussion of the theory of the decay of metastable phases (see, e.g., Refs. lT3T^) . but 
only give a brief summary. The KJMA theory in its simplest form describes the decay of a metastable phase to an 
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equilibrium phase, driven by a difference in free energy, from simple basic assumptions: negligibly small nuclei form 
at random and subsequently grow isotropically to larger droplets. An extended volume is defined as the sum of the 
volumes of all droplets, which grow and overlap freely, even upon collision. Since the nucleation events are completely 
random, they can also occur in the already transformed volume. The true transformed volume is the space covered 
by at least one droplet. Due to the randomness of the process, the relation between the extended volume fraction 
</? G xt and the (true) transformed volume fraction ip can be shown to bei^ 

l-ip = exp(-ip cxt ). (1) 

Assuming a constant radial growth velocity v and a constant nucleation rate / (corresponding to continuous nucle- 
ation), the extended volume is, in two dimensions, 

<^cxt = — g — , (2) 

at time t after the beginning of the phase transformation. Continuous nucleation is often associated with homogeneous 
nucleation due to thermal fluctuations. Another case, initial nucleation, where all the nuclei are present at the very 
beginning of the transformation, is often associated with heterogeneous nucleation at defects. In the case of continuous 
nucleation, thermal fluctuations that lead to the formation of small droplets obey metastablc thermodynamics, and 
knowledge of the size dependence of the free energy F of compact small droplets is needed. In a mesoscopic description, 
which should hold sufficiently below the critical temperature, this is obtained as a sum of a surface and a volume 
term. In two dimensions, this becomes 

F(r) = 2irra -(ji- Ho)A0nr 2 , (3) 

for a cluster with effective radius r and surface tension <jq along a primitive lattice vector at equilibrium (whose 
zero-field Ising analog is known exactly^). The difference in coverage AO between the two degenerate equilibrium 
phases corresponds to the difference of the spontaneous magnetization of the two degenerate phases at zero field in the 
Ising model»22, The use of the zero-field values is believed to be justified, and a recent publication^ supports the point 
of view that the surface tension is not strongly field dependent (though that study addresses a lower temperature). 
Conversely, the use of the macroscopic surface tension for microscopic structures is questionable, and more accurate 
expressions are a matter of research^SiSiS. and are sought in this paper as well. The free energy has a maximum 
at the critical droplet size (volume) s*. The existence of this maximum is the reason why the phase transformation 
is activated. On average, when smaller than s*, the droplets tend to dissolve, and when greater than s*, they tend to 
grow. The critical droplet size is a statistical quantity and is less sharply defined the flatter the shape of the maximum 
of the free energy. 

The velocity and the ratio of nucleation and growth are influenced by lateral mass transport. In ultra-high vacuum 
(UHV) and under the impact of a constant flux of the adsorbate, desorption may often be neglected (irreversible 
growth). In electrochemistry (EC), where the control parameter is the electrochemical potential, quasi-equilibrium is 
established for the subcritical fluctuations, so that desorption has to be included in a model description (reversible 
growth). As a consequence, in UHV, surface diffusion is generally the only means of lateral mass transport, while in 
EC, it can as well occur through the solution by independent desorption and adsorption events at different locations. 
In surface science, it has long been established that diffusion introduces correlations between clusters, since it reduces 
the nucleation probability in the vicinity of existing islands^2& Recent publications have addressed this violation of 
the basic assumptions of KJMA kinetics^Z^ Since the number of nucleation events is effectively diminished, phase 
transformation is slowed down. However, it has not been considered that the enhanced mass transport can increase 
the nucleation and growth rates, resulting in an acceleration of phase transformation. Thus, truly microscopic model 
calculations are desirable, and we provide them in this paper. This microscopic approach distinguishes our work from 
most other studies that aim at generalizing the KJMA theory. 

The rest of this paper is organized as follows: We introduce the model and describe our model calculations in Sec.lTTl 
Results are presented in Sec. IIIII First, we discuss the phase-transformation kinetics (Sec. IIII Afl and how they are 
affected by diffusion. Second, we discuss the effects of diffusion on the size distributions of supercritical (Sec. IIII B|) 
and subcritical (Sec. 1111 C[) clusters. For the latter, we give a theoretical explanation in terms of the free energy of 
lattice animals. A summary and conclusions follow in Sec. IIVI 
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II. MODEL AND ALGORITHM 
A. Model 

We describe the electrochemical deposition of ions onto a metal surface in terms of a square lattice-gas model with 
attractive nearest-neighbor interactions. The grand-canonical Hamiltonian of the system reads 

W = Ci(t) Ci (t)- a* X) <=*(*). (4) 

where <p is the interaction constant and [i the adsorbing ions' electrochemical potential in the solution, which represents 
an infinite reservoir. The q are occupation variables at site i with the values 




if empty 

1 if occupied. 



The first sum in Eq. Q runs over all pairs of nearest-neighbor sites, and the second runs over all lattice sites. Water 
and counter ions are not explicitly included. The Hamiltonian is equivalent to the Ising model, up to an additive 
factorM^ 

We perform model simulations for an L x L — 256 x 256 square lattice, applying periodic boundary conditions, 
for the values of the model parameters <f> = 4, and — Mo i = 0.4j2i where fio is the electrochemical potential at 
coexistence, and /i > no favors adsorption. The temperature T equals 0.8 T c w 0.453873 1.815348 J^S with T c 
the critical temperature of the Ising model for the given parameters^ Energy and temperature units are chosen such 
that Boltzmann's constant = 1. 

B. Algorithm 

To investigate the kinetics of the phase transformation, we perform kinetic Monte Carlo simulations )*■ where the 
transition probability for an elementary step is given by the Glauber acceptance probability 33 



1 + exp [(£/ - E^/foT] ' y ' 

with Ef and Ei the configurational energies in the final and the initial configurations, respectively. We choose 
this dynamic for the sake of consistency with earlier work^ Other dynamics that are commonly used in studies of 
electrochemical adsorption, such as the one-step dynamio 3 ^^ and the transition dynamic algorithm2S* 3 i also contain 
a local transition barrier and are co mmon ly known as Arrhenius dynamics. For further discussion and applications 
of such dynamics, see, e.g., Refs. Hll6l l38l39l 

Two kinds of elementary steps are considered: an adsorption/desorption step, which corresponds to a non-conserved 
order parameter, and a diffusion step to a nearest-neighbor position, which conserves the order parameter A new 
configuration is chosen by randomly picking a lattice site and proposing a trial step: either, with probability 1), 
an adsorption/desorption step or, with probability R/(R+1), a diffusion step to a randomly chosen nearest-neighbor 
position. Thus, the ratio of the diffusion attempt probability to the adsorption/desorption attempt probability equals 
R. Explicitly, we propose adsorption when we choose an empty site for an adsorption/desorption step, and desorption 
when we choose an occupied site - there is no second-layer growth. When we choose an empty site for a diffusion step, 
the initial configuration is maintained - vacancies do not diffuse. Only after attempting an adsorption/desorption step 
is the simulation clock updated by one Monte Carlo step (mcs); hence, a higher diffusion rate R does not slow down 
the time scale for adsorption/desorption, which is the physically relevant time scale for the phase transformation. 
After L 2 mcs, the simulation time thus has increased by one Monte Carlo step per site (MCSS). Hence, one MCSS 
means in our simulations one adsorption/desorption attempt per site. 

We start our simulations with an empty lattice, set the electrochemical potential to (i — /iq + 0.4, and let the 
system relax toward the stable configuration. We sample the time development of the coverage 9 (the appropriate 
order parameter of the system, equivalent to ip in KJMA theory) and, at certain coverages, evaluate the cluster-size 
distribution, using the Hoshen-Kopelman algorithm^ We perform model calculations for various values of R from 
to 1000. Every calculation is the average over 100 independent simulation runs. 

One of our goals is the understanding of the electrochemical deposition and subsequent dissolution of metal ad- 
sorbates after short pulses of the electrode potential. To that end, in addition to the constant-potential simulations 
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FIG. 1: Snapshots at different stages of the phase transformation during the potential-pulse simulations. R — 0. Coverages 
(from upper left to lower right): 8 — 0.05, 8 = 0.25 and 8 — 0.5 during the adsorption step, 8 — 0.25 and 8 = 0.05 during the 
desorption step. 



described in the previous paragraph, we perform potential-pulse simulations for values of R from to 10 000 starting 
from a metastable low-coverage configuration (result of an equilibration at /i = fiQ — 0.4), set fj, = fiQ + 0.4 and wait 
until a certain coverage S is reached. Then, we switch back to \i = fiQ — 0.4, at which the adsorbate desorbs again, 
and let the system equilibrate. 



III. RESULTS 



In this section, we present the results of our model simulations for the constant-potential and the potential-pulse 
cases. Snapshots of the adsorbate morphology at different time steps (Fig. show the nucleation, growth, and 
coalescence of clusters of quite irregular shapes (a more detailed discussion of this figure follows in Sec. lIIIBl"|) . It has 
been established, however, that the classical KJMA theory gives a good statistical description of phase transformation 
with microscopic dynamics, at least at its earlier stages and in the absence of diffusion, even though it assumes the 
growth of circular droplets^ Despite that the KJMA theory does not make assumptions about the nucleation and 
growth mechanisms, it is not clear how far it will hold when adsorbate diffusion is included in the kinetics of the 
lattice-gas model. We shall address this question in the following. 
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(f/T) 3 

FIG. 2: Time development of the coverage during constant-potential simulations, (a) Linear plot vs time, (b) Logarithmic plot 
of fraction of uncovered area vs cube of normalized time. 

TABLE I: Metastable lifetime for different diffusion rates R: Tg from 8 = 0.5 in constant-potential simulations and r p from 
peak position in potential-pulse simulations. 

R re (MCSS) r p (MCSS) 



187 184 

1 185 184 
10 179 177 
100 160 160 
1000 132 130 
10 000 111 



A. Phase-transformation kinetics 

Monte Carlo results for the time development of the coverage at constant potential [Fig. Eta)] apparently show the 
shape familiar from KJMA theory^ but with an acceleration of the phase transformation due to diffusion. For a 
diffusion rate of R — 10 the effect is not very marked, but it becomes quite pronounced for higher R. In order to 
quantify the effect, we define the metastable lifetime r as the time when the coverage reaches 9 = 1/2. The results 
are shown in Table [I] The lifetime r decreases weakly from i? = 0toi?=10, and then more strongly for larger R. 

KJMA theory predicts a linear dependence of the logarithm of the area fraction of the untransformed phase ln(l — 9) 
on t n with the Avrami exponent n = 3 [Eqs. and ©]. In order to distinguish a change in the linear dependence from 
the acceleration of the phase transformation, we normalize the time scales to the metastable lifetime r. Figure Etb) 



6 



TABLE II: Kinetic parameters of the KJMA theory from fits of the phase transformation plots similar to Fig.[5fb) for different 
diffusion rates R. I nucleation rate, v linear growth rate, r metastable life time. 



R Iv 2 (MCSS 3 ) Iv 2 t s 

9.1 x 1(T 8 0.59 

1 9.5 x 10~ 8 0.59 
10 1.0 x 10~ 7 0.59 
100 1.5 x 10" 7 0.59 
1000 2.5 x 10" 7 0.56 



shows almost perfect data collapse up to R = 100. At larger R, the slope is initially steeper. This corresponds to a 
relatively higher increase of the transformation rate, which levels off at t ss r. To obtain quantitative information on 
the kinetics, wc perform for all values of R unweighted least-squares fits to straight lines of the corresponding data 
points with the time scale unnormalized. The slopes give the product Iv 2 . From Eqs. Q and @, the contributions 
of the nucleation and growth rates, / and v, to the acceleration of the phase transformation cannot be distinguished. 
For a separate estimate of / and v, one needs to know the variance of the order parameter jAS. which we did not sample 
in our simulations. 

The fits have to be restricted to a certain time interval. Around t = t, the plots show a slight kink, a result of surface 
tension effects around the percolation coverage, which disturb the independent growth of the droplets. Therefore, we 
choose 0.9 x r as the upper limit. At very short times the system relaxes quickly to the metastable state. Therefore, 
we choose the lower time limit by stepwise decreasing it until x 2 per degree of freedom starts to increase steeply. We 
here refrain from correcting for the background of the metastable phase, in contrast to Ref. ^[ 

The results for Iv 2 shown in Table [H] reflect the acceleration of the phase transformation. The product Iv 2 is 
inversely proportional to the cube of the metastable lifetime r, and as long as the dynamics of the phase transformation 
do not change qualitatively, the situation is tantamount to a mere rescaling of the time scale. Then, the film 
morphology at the same coverage should not depend on R. The product Iv 2 t 3 is constant for R < 100, but differs 
at R = 1000, indicating a qualitative change in the phase-transformation dynamics for large R. Accordingly, we find 
that the curve for the same R fits less well to a straight line, so that the fitting interval is much shorter. We will 
return to this point in the next paragraph. In Sec. HUB 21 we will show that this change coincides with a change of 
the film morphology. 

The adsorption step of the potential-pulse calculations is equivalent to the constant-potential simulations. After the 
potential switch, the adsorbate desorbs, resulting in peak-shaped coverage curves, as we demonstrate for a switching 
coverage of 0.5 in Fig. CJa). The peak position gives another estimate for the metastable lifetime r (Table [|J. It 
coincides with the former estimate. In the absence of diffusion, the desorption is faster than the adsorption - the 
desorption can proceed without previous nucleation by simple shrinking of the clusters. The asymmetry in the kinetics 
between adsorption and desorption decreases with increasing R and is very small at R = 10 000 [Fig.^b)]. At this 
diffusion rate, the decay of the metastable phase seems to be much better described by an Avrami exponent n = 2, 
consistent with instantaneous nucleation (Fig.^J). 

During the growth of the adsorbate phase, nucleation may be effectively inhibited in the capture zones around 
clusters, and the size of these zones depends on the diffusion rate424£ Concerning the apparent Avrami exponent 
n = 2 for R = 10 000, one may speculate that nucleation is initially very fast, but gets inhibited quickly on the whole 
surface, which is completely covered by the capture zones of existing clusters. The result is instantaneous, though 
homogeneous nucleation. The subsequent adsorption kinetics are governed solely by the growth rate, as is the case for 
the desorption kinetics after the potential step. This explains the almost symmetric coverage curves for R = 10 000 
[Fig. Eft)]- 

B. Cluster-size distributions 

We characterize the film morphology by means of the cluster-size distributions. We show the number densities, 
i.e. the total number of clusters of size s on the finite lattice, divided by the number of lattice sites. This equals the 
probability of finding on a site the center of a cluster of size s. The size s is the number of particles in the cluster. 
Due to the low number densities for large cluster sizes, we use a growing bin size for the histograms, resulting in an 
equal distribution of data points on a logarithmic scale. We show results for the potential-pulse simulations. 
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FIG. 3: Time development of the coverage during potential-pulse simulations with switching coverage 9 a — 0.5. (a) Linear, (b) 
Mirrored at r and normalized with r. Curves labeled (ads) during adsorption step, curves labeled (des) during desorption step. 

1. In the absence of diffusion 

In the absence of diffusion, we show the evolution of the cluster-size distributions at different stages of the course 
of the potential-pulse simulations as log-log plots in Fig. [SJ During the adsorption step, the maximum cluster size 
increases with time, as existing clusters grow to larger sizes. The number densities of the cluster sizes that had been 
present previously decrease. This indicates that with time fewer clusters nucleate and grow into a certain size interval 
than grow out of it. In the early stages of the desorption step, it can be clearly seen that, while the densities of the 
largest and of medium-sized and small clusters are decreasing, some intermediately large clusters (around s — 1000) 
still gain in population. This is also true in the presence of diffusion (not shown). Later, the distribution gets consumed 
from the large-cluster side. Hence, the distributions at equal coverage differ markedly between the adsorption and the 
desorption step. This difference is clearly reflected in the snapshots of the film morphology in Fig.^ when comparing 
equal coverages during adsorption and desorption. The latter show fewer, but larger clusters. These snapshots also 
demonstrate the effect of coalescence during adsorption on the morphology - individual structures that form larger 
aggregates at later times can be clearly identified. Conversely, during desorption, dissociation of smaller clusters from 
larger aggregates may occur. A complete theory of cluster-size distributions, which we do not aim to present here, 
must consider these complications. We conclude that, during adsorption, existing clusters grow to larger sizes and 
may undergo coalescence, while new and thus smaller clusters can only appear on the untransformed area, thus with 
reduced probability. In the early stages of desorption, ripening of the morphology occurs in parallel with a decrease 
of the coverage, while later, clusters of all sizes shrink at constant rates. 
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FIG. 4: Time development of the coverage during potential-pulse simulations up to r. Logarithmic plot of fraction of uncovered 
area vs (a) cube (b) square of normalized time. The thin straight lines are guides to the eye. 
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FIG. 5: Cluster-size distributions (number of clusters of size s per unit area) for potential-pulse simulation during the adsorption 
step (ads) and the desorption step (des). Switching coverage 6 a = 0.5, diffusion rate R = 0. 
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FIG. 6: Cluster-size distributions (number of clusters of size s per unit area) for potential-pulse simulations with switching 
coverage 8 S — 0.5. (a) 8 = 0.05 (during adsorption step), (b) 6 = 0.5 (potential switch), (c) 9 = 0.05 (during desorption step). 

2. In the presence of diffusion 

For the potential-pulse simulations in the presence of diffusion, results for different diffusion rates at 9 = 0.05 during 
adsorption, 9 = 0.5 (the switching coverage), and 9 = 0.05 during desorption are shown in Fig. as log-log plots. In 
the early stages of growth, we find that for fast diffusion rates (R > 1000), the maximum cluster sizes present are 
larger [Fig. [fjfa)]. At later stages, there is a marked reduction through fast diffusion of the densities of medium-sized 
clusters, together with an increase of the densities of large clusters [Fig. E{b)]. In accordance with this observation, 
snapshots of the phase transformation show a coarser morphology for R = 1000 than for R — (Fig. EJ. Near the 
end of the reverse step, the reduction of the amount of medium-sized clusters prevails, but much weaker, such that 
the morphologies for different R get more similar during desorption [Fig. EJc)]. For R < 100, the distributions are 
virtually indistinguishable from the case without diffusion. This indicates that the phase-transformation dynamics 
are similar and just proceed on a faster time scale (Sec. IIIIA"fr . 

We find the reduction of medium-sized clusters more likely to be a result of increased coalescence rates than of a 
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FIG. 7: Snapshots for diffusion rates 7? = (left) and R — 1000 (right). In both cases, 9 = 0.5. 



TABLE III: Reduction of the monomer number density n\ relative to the untransformed surface fraction. R diffusion rate, 
coverage. See Sec- EHnH 



R 



m(0) 



1-8 



»l(g) 
n 1 (0.O5) 



1000 



0.05 

0.25 

0.5 

0.05 

0.25 

0.5 



1.25 x 10~ 2 
9.10 x 10" 3 
5.53 x 10" 3 
1.24 x 10~ 2 
8.60 x 10" 3 
5.51 x 10~ 3 



1 

0.79 
0.53 
1 

0.79 
0.53 



1 

0.72 
0.44 
1 

0.69 
0.41 



change of the radius-dependent growth law of the clusters. While the former can be a result of a collective motion of 
clusters established by edge diffusion or of dissociation and diffusion of monomers and is well established in surface 
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the latter would imply an unphysical slowing-down of growth only at medium sizes. 



We have seen that diffusion affects the cluster-size distributions only at diffusion rates high enough that also the 
phase-transformation kinetics cross over from continuous to initial nucleation, R > 1000 (Sec. IIII*A*jl . Lower diffusion 
rates speed up the dynamics without changing them qualitatively. A rough estimation considering the number densities 
of monomers suggests that at these lower rates, diffusion is still not more frequent than adsorption/desorption, as the 
mere numerical values of R naively might suggest. 



C. Subcritical cluster distributions 



Unlike the supercritical droplets, whose size distribution is governed by the phase-transformation dynamics, the 
subcritical droplets can be considered to be quasi-equilibrium fluctuations. This allows a thermodynamic approach for 
a theoretical description of their size distribution. We first show simulation results of the size distributions. Second, 
we calculate their theoretical quasi-equilibrium distribution within a lattice-animal model. 



1. Simulation results 



We show the number densities, as explained in Sec. IIII Bl of the subcritical clusters for potential-pulse simulations 
with a switching coverage of 9 S = 0.5 in Fig. |SJ At all stages of the simulation, the size distribution of small clusters 
decreases monotonically with increasing size. The distribution develops with time. During the adsorption step, the 
number densities decrease. This decrease is somewhat greater than expected from the decrease of the untransformed 
surface fraction. Column 4 in Table ITTTl shows the fraction to which the free surface gets reduced compared to 9 = 0.05, 
and column 5 shows the corresponding reduction of the monomer density. The reduction of the monomer density is 
always greater than the reduction of the free surface, and it is enhanced by increasing diffusion rate. This indicates 
that in a capture zone close to an existing cluster nucleation is suppressed, since aggregation to the cluster is more 
likely. The size of the capture zones is increased by diffusion. During the desorption step, the number densities of 
small clusters are generally further diminished (with the exception of monomers and dimers, whose densities increase). 
Again, in the course of desorption, the distributions with and without diffusion approach each other. 
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FIG. 8: Subcritical cluster-size distributions (number of clusters of size s per unit area) for R = (circles) and i? = 1000 
(squares) at different stages of the potential-pulse simulations. Curves labeled (ads) during the adsorption step, curves labeled 
(des) during the desorption step. 

2. Lattice- animal model 

The cluster-size probability distribution gives the probability P{s) that a subcritical cluster is of size s, and it is 
normalized such that the sum of all probabilities equals 1. To calculate it, we need an expression for the free energy 
F(s) of a cluster of size s, since 



P(s) oc exp(~F(s)/k B T). 



(6) 



The theoretical prediction from droplet theory [Eq. 0], in two dimensions, using \fs~pK as an effective radius r and 
the exact zero-field Ising values for uq and AC, is shown in Fig. |5J The agreement is rather poor. Since the surface 
tension, which is a macroscopic quantity, may not be a good description for such small clusters, we seek a better 
statistical-mechanical expression for the cluster free energy. If we knew all possible configurations of a cluster of size s 
and their numbers of bonds to nearest-neighbor unoccupied sites (including holes) t, we could compute its restricted 
partition function exactly as 



Z(s) = cxp 



s (M ~ Mo) 



t 



w(s, t) exp 



tip 



2 kWT 



(7) 



More specifically, we need to know the number w(s,t) of degenerate clusters with same s and t. These numbers are 
very high unless s is very small. As a consequence, even with the help of computers, they are accessible only up to 
moderate values of s. Generating these configurations is known as counting lattice animals ^^JO^L^H In the lattice- 
animal literature (see, e.g., Ref. l5lf) . the perimeter, if considered, is usually defined as the number of unoccupied 
nearest-neighbor sites and not as the number of bonds to these sites, and only w(s) is known. For the calculation of 
w(s, t), we use the algorithm of Redner^ We show the results for animals up to s = 21 in Table Hvl The computation 
took about 65 hours of CPU time on an Intel Xeon-based workstation. The w(s) from our calculations are in accord 
with the known rcsults^SiSl From the partition function, the free energy is calculated as 



F(s) = -k B ThiZ(s), 



(8) 



so that 



P(s)=Z(s)/J2z(s), 



(9) 



summing up to the critical cluster size s* . This calculation, as the droplet theory, assumes non- interacting subcritical 
clusters which are in quasi-equilibrium throughout the course of the phase transformation, an assumption that should 
be better fulfilled at higher temperatures, as is the case in our simulations. The resulting cluster-size probability 
distribution is shown and compared to simulation results in Fig. |5| At early times and in the absence of diffusion, 



12 




10 



10 



G O lattice animals \ 

□ □ droplet theory ' ■ ^ 

MC, 9 = 0.05, R = ' s . 

MC, 9 = 0.75, R = 

MC, 9=0.75,fi=1000 

MC, 9=0.05 (des), R = 

_j i i i I i i i i 



5 10 
Cluster size s 



15 



FIG. 9: Subcritical cluster-size probability distributions (probability that a subcritical cluster is of size s): Monte-Carlo results 
and theoretical predictions from droplet theory and from counting lattice animals. Curve labeled (des) during the desorption 
step, others during the adsorption step. 
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FIG. 10: Free energy of small clusters according to droplet theory and from counting lattice animals. The critical cluster size 
is labeled s*. 



the agreement is reasonable (note the logarithmic scale). The trends that we observed for the number densities are 
reflected as an increasing reduction of the probabilities with increasing cluster size for later times, in the presence 
of diffusion, and during desorption, respectively. We can attribute these trends to increasing deviations from the 
quasi-equilibrium distribution - the larger subcritical clusters are consumed faster than they are produced, or during 
desorption, the small clusters are produced faster than they disappear (remember their increase in absolute density). 
In contrast, the predictions of droplet theory are way off. This is a consequence of the surface-tension contribution 
to the free energy of the clusters (Fig. I10|) ; surface fluctuations from the continuum description, which contribute 
to the partition function, are unphysical on the discrete lattice for small clusters. In the lattice-animal model, the 
maximum of the free energy is shifted to larger clusters compared to droplet theory, thus giving a critical cluster 
size of 18 instead of 12. In both cases, the free energy changes very little in the neighborhood of the critical cluster. 
An approach analogous to the one presented in this section was independently developed in Ref. |25| (with w(s, t) 
computed up to s — 17). 



IV. SUMMARY AND CONCLUSIONS 



We have performed kinetic Monte Carlo simulations of a square lattice-gas model using Glauber dynamics to 
investigate the decay of a metastable low coverage phase. The model, which contains adsorption/desorption and 
nearest-neighbor diffusion moves, is used to qualitatively describe the electrosorption and subsequent dissolution of 
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TABLE IV: Number w of configurations of lattice animals of size s and perimeter t. For the definition of the perimeter, see 
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metal ions onto a metal substrate in potential-pulse experiments. We have observed the phase-transformation kinetics 
and the size distributions of sub- and supercritical clusters. The results have been compared to the KJMA theory. 

Nearest-neighbor diffusion speeds up the phase transformation. At low and moderate diffusion rates, the kinetics are 
well described by the KJMA theory for continuous nucleation. At high diffusion rates, we find a crossover to kinetics 
that resemble initial nucleation. In the former case, the cluster-size distributions are virtually indistinguishable from 
each other, but in the latter case, the densities of medium-sized clusters get reduced in favor of large clusters. At the 
same coverage, the morphology is coarser during the desorption step than during the adsorption step. The crossover 
to initial nucleation is likely to be a consequence of the diffusion-rate dependence of the size of the capture zones of 
clusters. When the diffusion rate increases, the capture zones get larger, and for very fast diffusion rates, they might 
become space-filling already in early stages of the phase transformation. 

We have calculated the free energies of subcritical clusters from their exact restricted partition functions, which we 
obtain from counting lattice animals. This enables us to calculate a theoretical prediction of the subcritical cluster- 
size probability distribution. The prediction is in reasonable agreement with simulation results, apart from some well 
understandable deviations from the quasi-equilibrium distributions. In contrast, classical droplet theory fails, since 
the macroscopic surface tension is a poor description for the discrete microscopic clusters. 

Berthier et all— use a similar way as we for the calculation of subcritical cluster distributions. To obtain the 
number density of clusters of size s, they consider the configurational entropy of the ground and the first two excited 
states. These are the clusters with the three lowest perimeters t. Compared to our partition function, which considers 
all possible configurations (up to the highest excited state), their summation is incomplete for s > 7. Moreover, 
many of their weights for the excited states are too low. Since they do not specify the method of counting the 
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accessible configurations, the origin of the discrepancy is unclear, but it seems that they miss some configurations. 
As a consequence, the free energies obtained by their method are erroneously non-monotonic and too high. 

The experimentally observed post-deposition dynamics after a potential pulsep? which result in the continued 
growth of large clusters while small clusters already tend to dissolve, are clearly reproduced by our simulations. 
Though it is not necessary to explicitly include surface diffusion under control of the electrochemical potential to 
establish lateral mass transport, we find diffusion more likely to be the dominant mechanism under the experimental 
conditions. However, in our simulations, coverages are much higher than in the experiment. While time scales of 
adsorption and desorption in the experiment are rather different, in our model they are necessarily much more similar, 
since adsorption and desorption proceed with very similar mechanisms. 

Our model calculations are a first step toward including diffusion into the microscopic dynamics of a kinetic lattice- 
gas (or Ising) model and to investigate its influence on the applicability of the KJMA theory. Still, calculations over 
a larger range of fields and a more stringent analysis are required^ To support our hypotheses about the interplay 
between diffusion, nucleation, and spatial correlation, more microscopic information, like on nucleation events and the 
fate of individual clusters, should be obtained from the simulations. For a more meaningful description of experiments, 
a closer reproduction of experimental conditions is desirable. 
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